Short-range interactions between strongly nonlocal spatial solitons 
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A novel phenomenon is discovered that the short-range interaction between strongly nonlocal 
spatial solitons depends sinusoidally on their phase diilerence. The two neighbouring solitons at close 
proximate can be inter-trapped via the strong nonlocality, and propagate together as a whole. The 
trajectory of the propagation is a straight line with its slope controlled by the phase difference. The 
experimental results carried out in nematic liquid crystals agree quantitatively with the prediction. 
Our study suggests that the phenomenon to steer optical beams by controlling the phase difference 
could be used in all-optical information processing. 
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Solitons are a common phenomenon appearing in many 
physical fields, while the interactions of the solitons 
have great potential for much wild applications, [ij] The 
strongly nonlocal spatial soliton, which is also known 
the accessible soliton [2"| and is the self-trapped optical 
beam from the balance between diffraction and nonlin- 
earity propagating in nonlocal nonlinear media under the 
condition of the strong nonlocality [3], has some significa- 
tive properties and has attracted more attentions in the 
last decade Several strongly nonlocal, referred 

also as highly nonlocal in some literatures (for example, 
Refs.i, 11,11), media have been found in experiments, such 
as the nematic liquid crystal [3,[1,[^01, the lead glass @, 
and the thermal nonlinear liquid [9"], and the nonlinear 
gas of ions[l(|. 

Nonlocality of nonlinear response may drastically mod- 
ify the properties of the solitons, specially their in- 
teractions. In a strongly nonlocal case, it has been 
shown theoretically and experimentally that the attrac- 
tion can occur between the bright solitons with any 
phase difference [il. [ill [l^. ITs^ . the coherent or incoherent 
solitons nil, or the dark solitons [^, [3]. In the con- 
trary, however, the interaction can be attractive only for 
two in-phase local bright solitons^. Both of the long- 
range interaction [l^ and the short-range interaction be- 
tween the solitons can happen in strongly nonlocal non- 
linearity, while only the short-range interaction can occur 
in local nonlinearity because the force between the local 
solitons decreases exponentially with the separation be- 
tween them [18 1 . 

The fact well-known so far was that3,[il[i3,[illi3 

the interactions between the strongly nonlocal bright soli- 
tons are independent of their phase difference, in this 
letter, we differentiate the patterns of the short-range 
interaction and the long-range interaction between the 
strongly nonlocal bright solitons. We show theoretically 
and experimentally that the short-range interaction of 
the two strongly nonlocal bright solitons is sinusoidally 
dependent on their phase difference. 

Let us consider a (l-|-2)-D model of an optical field that 
polarizes linearly with an envelope A and propagates in 



z-direction in the medium with nonlocal nonlinearity 
2ikd^A + \J\A + 2fc2 A = 0, 



(1) 



where V'^ = 9^ + 9,^ , k and ng are the wave- vector and 
linear refractive index of the medium. The nonlinear per- 
turbation of refraction index An(a;, y, z) can be generally 
expressed as 



An = n2 R{x-x',y-y')\A{x',y',z)\^dx'dy', (2) 

where R{x, y) is the real nonlinear response func- 
tion of the medium. The normalized condition, 
R{x, y)dxdy — 1, is chosen physically to make the 
nonlinear index n2 have the same dimensions as that 
in a local Kerr-medium. If R{x,y) is a delta function. 
An = n2|/lp and Eq.(IT]) becomes the well-known non- 
linear schrodinger equations (NLSE) for the local Kerr- 
medium. The Eq.(IT|) and Eq.(I2), so-called the nonlocal 
nonlinear schrodinger equation (NNLSE) , can model the 
beam propagation in most of nonlocally nonlinear media 
discovered in experiments so far. 

For the NNLSE, several well-known invariant integrals 
are important for the theoretical analysis [l^ The 
first one is the power integral, P = \A{x,y)\'^dxdy, 
which results from the energy conservation of an optical 
beam during the propagation in a lossless medium. The 
second is the momentum 12 , loj 



M 



A*{-'iV±)Adxdy. 



(3) 



The momentum governs the movement of the mass center 
of the two optical beams, i.e. 



dz 



M 
"P ' 



(4) 
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where the mass center r 

(1/P) / f^A(x,y,z)\'^dxdy, and f — xe^ + ycy. The 
trajectory of mass center is a straight line with its slope 
determined by Eq. 



rdz) 
f ycy. 
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FIG. 1: The sketchs of (a) the two injective solitons, (b) 
the hquid crystal cell, and (c) the experimental setup. NA, 
neutral attenuator; BS, beam splitters; M, plate mirror; PP, 
parallel-face plate for adjusting the phase difference; O, lOx 
microscope objective; LC, liquid crystal cell; MS, microscope; 
F, laser-line filter; BP, beam profiler. 
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FIG. 2: The dependence of the slope on the distance h (a) 
and phase difference 7 (b) for two parallel-injected solitons. 
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FIG. 3: The contour graph of the numerical propagation of 
Eqs. IT} and ((2} for the two parallel-injected solitons. The 
local case is shown in the first (from the left to the right) col- 
umn, and the two different nonlocal cases with an exponential- 
decay response function given by Eq. (O are in the second 
and third columns, respectively. The red solid lines show the 
movement of the mass center of the two solitons. The re- 
sults from the Snyder- Mitchell mode! 2] with the same initial 
condition are also given in the fourth column for comparison. 
The phase difference between two solitons are 0, 7r/2, vr, and 
37r/2, respectively (from the top to the bottom). 



Suppose the two simultaneously incident Gaussian soli- 
tons are coplanar in x — z plane, with a width wq, a 
phase difference 7 and a separation d{^— 2h), as shown in 
Fig.IHa), i.e. 



A{x,y,0) = Aoexp 

+ Aqc^'^ exp 



(x + hf + 



-f ik{x + h) tan a 

{x - hY + 



2wl 



2wl 



ik{x — h) tanc( 



where the amplitude Aq is large enough to make the two 
beams propagate in soliton states [^, [llj. For the 
input condition ([5]), the slope can be obtained 



tan 


{h/wo) exp 


— (h/wo)'^ — (tana/9)^ 


sin 7 


e 


1 -|- exp 




- (h/wo)'^ — (tana/O)^ 


cos 7 



(6) 

and tanPy = 0, where tan f3x and tan/3j, are the slops in 
X — z and y — z planes respectively, and = 1/ kwQ is 
the far-fields divergence angle of a Gaussian beam. 

The slope of the line for the trajectory of mass center 
is greatly dependent on the separation d and the phase 



difference 7, as shown in Fig. [21 where we take a = 0, 
that is, two solitons are parallel injected into the medium. 
Figure [2] shows that tan /3x = only when 7 = or 
TT for h/wQ < 2, and tanfix goes toward to zero when 
h/wo > 2. tan has significant value when h is about 
or smaller than the beam width wq. In other words, 
when the soliton separation d is approximate or larger 
than four times the soliton width wq, the optical fields of 
) the two solitons do not overlap so that tan jSx decreases to 
zero, otherwise the two solitons are of an effective overlap 
and tan Px has a non-zero value changeable by the phase 
difference. 

It is important to emphasize that the above analytical 
result about the movement of the mass center [Eq. ([6])] is 
universal, independent of the form of nonlinear response 
function R. This means that no matter what the material 
is and how the degree of nonlocality is, the movement of 
the mass center is the same for the initial condition ([5]). 

Although the momentum gives the movement of the 
mass center, it is difficult to obtain the analytical so- 
lution of the beam propagation for the initial condition 
We carry out the numerical simulation for local, 
weakly nonlocal, and strongly nonlocal propagations, re- 
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spectively, and only the (1+1)-D case of Eqs. ([T]) and 
([2]) is simulated for the sake of simplicity and also with- 
out the loss of generality. The (1+1)-D model makes 
it possible to compare the propagations in the nonlocal 
nonlinearity and in the local nonlinearity, and also pro- 
vides an exact enough description to the (l-|-2)-D copla- 
nar propagation. The results are shown in Fig. [3l The 
nematic liquid crystal (NLC) is taken into consideration 
here for the consistency with our experiment discussed 
later in this letter. The response function of the NLC is 
the exponential-decay function 



R{x) = {l/2wm)exp{-\x\/ww), 



(7) 



for the (1+1)-D case[lj| 
Bessel function, R{x, y) 



and the zero-th order modified 
= (1/2™^ 



for the (l-f2)-D cylindrical symmetrical case[J, He 
where Wm is the bias-voltage controllable characteristic 
length [l^ of the response function R for the NLC. The 
ratio Wm/wQ indicates the degree of nonlocalityQ, which 
are chosen in simulation to be 0.47 and 10 for weakly non- 
local and strongly nonlocal cases in Fig. [31 respectively. 

Some interesting consequences can be obtained from 
Fig. [3l For the local system described by the NLSE, the 
two solitons attract each other only for the in-phase case 
(7 = 0) and repel each other for other occasions [ij. There 
also exists the power transfer between the solitons when 
7 is not equal and tt, as mentioned in Ref. [ij. The 
force between the in-phase solitons is always attractive, 
independent of the degree of nonlocality, as shown in the 
first row; while the repelling force between the solitons 
for the other phase difference cases becomes weaker, as 
the degree of nonlocality increases. As a result, the two 
solitons with an arbitrary phase difference can get attrac- 
tive when nonlocality becomes strong enough. For all of 
those propagations, however, the movement of the mass 
center obeys the same regulation, a straight-line trajec- 
tory with its slope given by Eq. ([6]). It is clear that the 
strong enough nonlocality can make the two spatial soli- 
tons trap each other and propagate together as a whole, 
going along the mass center trajectory, the straight line, 
with the slope steered by their initial phase difference. 
Interestingly, a similar phenomenon was mentioned (20j 
in the intermediate process |21|] for dealing with incoher- 
ent solitons in "fast" nonlocal nonlinear media. 

It can be observed that the maximum value of the 
tilting angles occurs when 7 approximates tt for small 
distance h. For each h, two values "fmax between 
7r/2 and 37r/2, which is determined by cos7maa; = 
± exp(— /i^/wq), makes f3x reach the extrema (a maxi- 
mum and a minimum, respectively), ta,n[{Px) e,j;t]/Q — 
±{h/wo) exp(— /i^/wo)/^/! — exp(— 2/i2/w2). The maxi- 
mum is for steering right and the minimum for left. The 
smaller h, the larger tilt angle Px- The largest tilting an- 
gle is e/v^. It means the steering angle of whole beam 
are significance only for thin beams. 
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FIG. 4: (a) the phase difference 'ymax for the maximum slope 
angle vs the distance h, and (b) the maximum of tilt angle 
tan[(/32,)ea^t] vs the distance h. Both are for q = 0. 



Togetherness of Figs.[2]and[3]can result in one more sig- 
nificant outcome. When d < Awq, the two solitons have 
a non-zero overlap and the slope tan (3x is no-zero also. 
In this case, two solitons can be inter-trapped via the 
strong nonlocality, and propagate together as a whole, 
going along the diagonal line of the trajectory for their 
mass center. This is the short-range interaction between 
the strongly nonlocal solitons, which is phase-sensitive 
(controllable by their phase difference). When d > 4wo, 
on the other hand, the two solitons never overlap and the 
slope of the trajectory for their mass center tends toward 
zero. In this case, two strongly nonlocal solitons undergo 
periodic collisions in the coplanar propagation [3. fill lrn| . 
and spiral about one another if they are initially skew 
to each other Both processes have nothing to do 

with their relative phase, as predicted first b y S nyder and 
MitcheUfi] and verified experimentally [Uj, |17|. This is 
the long-range interaction between the strongly nonlocal 
solitons. As the soliton separation d increases from less 
to larger than 4wo, the interaction will gradually tran- 
sit from the short-range pattern to the long-range one in 
the strongly nonlocal nonlinearity, and vice versa. Only 
the short-range interaction, however, exists in the local 
nonlinearity [18]. 

Two more points are given in the end of the theoret- 
ical part. It is worth to note that the Snyder-Mitchell 
model[2|| can not give right prediction about the short- 
range interaction between the strongly nonlocal solitons, 
as shown in the fourth column of Fig. [3l Such a phase- 
controllable short-range interaction might have its appli- 
cation in all-optical switching and routing. 

In order to verify our prediction, we carried out the 
experiment in the NLC. 

The configuration of the NLC cell is the same as in 
the previous works 0, 0, [l3|) as shown in Fig. [T](b). The 
optical field polarized in j/-axis with envelope A propa- 
gates in z-direction. An external low-frequency electric 
field EfiF is applied in y-direction to control the initial 
tiff angle of the NLC. 

The experimental setup is illustrated in Fig. [Ijc). 
The laser beam from the laser is split into two beams, 
then they are combined together with a small separa- 
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tion through the other beam-spUtter and launched into 
a 80/Ltm-thick NLC cell by a 10 x microscope objective. 
The beam width at the focus Wq, the separation d, and 
relative angle 2a between the two beams are measured 
by an edged-scanning beam profiler when the NLC cell 
is removed. The cell is filled with the NLC TEB30A 
(from SLICHEM China Ltd.), whose n|| = 1.6924, nj_ = 
1.5221, K « 10~"7V, e°P = 0.5474, and e^^ = 9.4. The 
bias voltage on the cell is set to 1.4V, and then a pretilt 
angle is nearly 7r/4 in order to obtain strong enough non- 
locality and the lowest critical power of sohtonsj^]. The 
launched power for each beam is fixed to 6mW, and two 
spatial solitons are obtained for such high enough exci- 
tations. The parameters for the beams inside the NLC 
are calculated from the measurement without NLC cell, 
i.e. Wo — 2.2/im, d — 2.25/im, tan(2Q!) = 0.0076, and the 
divergence angle 9 = 0.0231. 

The phase difference between the two beams (solitons) 
is adjusted by the rotation of a l.Smm-thick parallel- 
face plate, and measured through the interference pattern 
by the beam profiler located on the other branch after 
the second beam-splitter. First we find the position of 
the plate while the phase difference is adjusted to (in 
phase) , then we rotate the plate in small steps to increase 
the phase difference 7. 

We record the soliton trajectories for the different sit- 
uations by the CCD camera, as shown in FigO In Fig. [5] 
(a) and (b), each of two solitons is alone launched into 
the NLC respectively, and their trajectories are straight 
and horizontal. When two solitons are injected simulta- 
neously into the NLC, they will propagate as a whole, and 
tilt (c) up or (d) down (actually in a;-direction) . Since 
the separation is so small that two solitons cannot be dis- 
tinguished by the microscope in our experiment, we will 
see a whole beam, as a bound state, steered by the phase 
difference 7. 

In order to compare quantitatively our experimental 
observation with our theoretical prediction, we give the 
variation of the tilting angle with the phase difference 7 in 
Fig. [6] For each 7, we take five photos of the beam tilting 
angles to minimize the jitter of the tilting angle resulting 
from that of the laser source and the phase difference. 
We can see that the experiment points locate around the 
theoretical prediction with a relative small random er- 
ror. The error may mainly come from the slight jitter 
of the phase difference 7. Except those random error, 
we can say that the experiment results consist with the 
theoretical prediction very well. The maximum tilt angle 
observed in experiment is about 1.2°, approximate 0.6O. 

In conclusion, we predict that the short-range interac- 
tion between the strongly nonlocal solitons depends dras- 
tically on their phase difference. The result is universal, 
independent of the different nonlocal nonlinear media. 
The experiment carried out in the nematic liquid crystal 
agrees quantitatively with the prediction. The Snyder- 
Mitchell model can only give the right prediction for the 




FIG. 5: Photos of the beam trajectories for the single soli- 
ton [(a) and (b)] and the two solitons injected together [(c) 
and (d)] propagating in the NLC cell. The phase difference 
between the two solitons for (c) and (d) are about it/ 2 and 
3tt/2, respectively. 
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FIG. 6: The tilting angle of two beams vs. the phase differ- 
ence between them. Square points; experiment results, solid 
curve: the theoretical fitting curve from Eq. ((6}. 

long-range interaction between the strongly nonlocal soli- 
tons. 

Unlike their local counterpart, the strongly nonlocal 
solitons can exhibit both the short-range interaction and 
the long-range interaction. The two kinds of interactions, 
however, have different patterns, the former is phase- 
sensitive and the latter is not. Therefore, each of the two 
phenomena would provide a means of controlling light 
with light, and be thus potentially useful in developing 
all-optical signal processing devices. They are, of course, 
supposed to be applied in different situations. 
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